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We propose a new realization of softly broken supersymmetric theories as theories defined on stochastic 
superspace. At the classical level, the supersymmetry breaking is parameterized in terms of a single (in general 
complex) mass parameter, ^, describing the stochasticity of the Grassmannian superspace coordinates. In the 
context of the standard model with stochastic supersymmetry, the structure of the soft breaking terms has various 
characteristic features that can be tested in LHC experiments. Namely, at the classical level, the parameter, 
the universal soft trilinear coupling Ao, the universal gaugino mass mi/2 and the universal scalar mass mo are 
given solely in terms of ^; there are no other arbitrary parameters. The relations are — , Aq = 2^*, 
m-i/2 ~ ICI/2 and mo — 0. At the quantum level, these relations hold at a certain scale A which is a second 
free parameter. The soft scalar masses, zero at tree-level, are induced radiatively through the renormalization 
group equations at one-loop. With this pattern of soft breaking terms, large supersymmetric contributions to 
FCNC processes are avoided. As a concrete illustration of the proposed formalism, we consider a minimal 
model, which is just the constrained MSSM with the stochastic superspace relations amongst the soft-breaking 
parameters imposed at the scale A. We show that this theory is phenomenologically viable for a certain region 
in the A) parameter space. Some sensible extensions of the minimal model are then briefly discussed. 



I. INTRODUCTION 

Supersymmetry is the unique non-trivial extension of the relativistic Poincare symmetry of spacetime. It leads to field theories 
with improved ultraviolet behavior. Because of this, supersymmetry might have direct relevance to particle physics through its 
stabilizing of the electroweak scale under radiative corrections. If this is indeed the case, then the non-observation of supersym- 
metric particles means that supersymmetry manifests at low energies in softly broken form: only those supersymmetry breaking 
operators are allowed which do not introduce new types of divergences (i.e. higher than logarithmic). 

Among renormalizable operators there are only three types of soft supersymmetry breaking terms: the mass term for the 
lowest scalar component of a chiral superfield, the mass term for the lowest fermionic component of a vector superfield and the 
scalar trilinear and bilinear couplings. In the supersymmetric version of the standard model, the generic supersymmetry breaking 
sector contains a large number of new unknown parameters. Consequently, the predictive power of the theory is compromised. 
Furthermore, there are in general unacceptably large contributions to flavor changing neutral current processes due to the absence 
of a GIM-like mechanism in the sector of scalar partners (squarks) of ordinary quarks. To avoid this problem the squark soft 
breaking masses must be degenerate with high accuracy, unlike the observed hierarchical structure of quark-lepton masses. To 
achieve universality in soft breaking masses is not an easy task. This typically requkes rather complicated hidden and messenger 
sectors of supersymmetry breaking and, in many cases, special flavor symmetries also need to be postulated. 

In this paper we suggest a conceptually different view of softly broken supersymmetric theories. Namely, we define a field 
theory on superspace where the Grassmannian coordinates are stochastic. With a suitably chosen probability distribution, to 
be defined below, the soft supersymmetry breaking parameters emerge upon averaging over the fluctuating superspace coordi- 
nates. At the classical level, the supersymmetry breaking is parameterized through a single mass scale, ^, that describes the 
stochasticity of the Grassmannian coordinates, and a very specific pattern of soft breaking terms emerges. Namely, at tree-level, 
the Bf^ parameter, the universal soft trilinear coupling Aq, and the universal gaugino mass mi/2 are given solely in terms of 
the stochasticity parameter ^, with no other arbitrary parameters involved. At the quantum level, these relations are imposed 
at a certain scale A, which becomes a second free parameter. The soft scalar masses, zero at tree-level, are induced radiatively 
through the renormalization group equations at one-loop. The resulting pattern of supersymmetry breaking is highly desirable 
phenomenologically since it is capable of resolving the above-mentioned supersymmetric flavor problem and it is extremely 
predictive. 

Most phenomenological studies of supersymmetric extensions of the standard model are done using explicit soft supersym- 
metry breaking. It is understood that one eventually wants to derive the soft breaking parameters from a fundamental theory 
of spontaneous or dynamical supersymmetry breaking, and to have that theory simultaneously explain the special pattern of 
soft breaking terms required for agreement with experimental constraints (degenerate squark masses, and so on). The idea of 
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stochastic superspace coordinates lies conceptually somewhere in-between a full fundamental theory and pure phenomenologi- 
cal parameterization. It is much more specific than the latter, but it is not intended to be a fundamental theory either. Our main 
point is that a certain line through the multi-dimensional parameter space of soft supersymmetry breaking turns out to be cor- 
related with a non-trivial observation: it corresponds to the hypothesis that superspace coordinates fluctuate and thereby break 
supersymmetry. Furthermore, a section of this line produces phenomenologically realistic physics, and being so constrained it 
is ripe for experimental confirmation or falsification. If the LHC were to discover a supersymmetric particle spectrum in accord 
with the stochastic superspace relations, that would motivate a new perspective on the desired fundamental theory. That theory 
would have to produce fluctuating superspace coordinates at the effective low-energy level. We shall not speculate here as to 
what such a theory could be hke. We begin at the beginning, namely phenomenology, and await the experimental results to come 
in the near future. 

The rest of the paper is organized as follows. In the next section we define stochastic superspace and as a warm up exercise 
compute soft breaking parameters in the simple interacting Wess-Zumino model. In the section 3 we apply our formalism to the 
case of the minimal supersymmetric standard model (MSSM). We compute the low-energy soft breaking parameters using an 
approximate analytical solution to the renormalization group (RG) equations, and establish that this model is phenomenologically 
acceptable in a certain region of A) parameter space. We briefly mention possible extensions and then conclude. 

II. STOCHASTIC SUPERSPACE 

We consider N=l 8-dimensional superspace, with the Grassmannian coordinates 9 and 9 being stochastic variables. One 
starts with the ordinary N=l supersymmetric Lagrangian density expressed in terms of superfields. The Lagrangian density 
in ordinary spacetime is then defined as an average of the superspace Lagrangian density over the different realizations of the 
stochastic Grassmannian coordinates described by some probability distribution, V{9, 9). This probability distribution can be 
expanded as': 

V{9, 9) = A + r *a + ^aS" + 6i"6i„B + 9j"C 
+9°'a^' Pv. + 9"9j^K^ + 9j'^9°'J:^ + 9"9j'^hD , (2.1) 

where, A, B, C, D and are c-numbers, while 'i', S, A and E are Grassmann numbers. We demand that the probability 
measure ( 12. Il l satisfy the following conditions: 

1 . The normaUzation condition: 

' d^9S9V{9,9) = 1, =^ D^l. (2.2) 



2. That all Lorentz non-scalar moments vanish to ensure Lorentz invariance: 

{9) = {§) = {99) = (9^9) = (99^) =0, =^* = S = V^=A = S = 0. (2.3) 
The non- vanishing moments are: 



^B^C* = -^, A=^. (2.4) 

where ^ is a complex parameter of mass dimension.^ 
Under these conditions we have the Hermitian probability measure, 

r{9, 9)\i? = 'PiO. 9)^1 + C'{00) + i{99) + \i\''{99){99) , (2.5) 



' We use conventions adopted in Ref. fl']. 

^ In tlieories with exact U(l)fl symmetry the phase of g has no physical meaning and can be rotated away. Below in numerical estimations we simply take 5 to 
be real. 
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for the stochastic Grassmann variables. We shall use the dimensionless probability distribution V in what follows. The prob- 
ability distribution in ( 12.51 ) can be considered as the analog of the more familiar white noise distribution for bosonic stochastic 
variables. It is convenient to introduce P{6) chiral and P{6) = P+ anti-chiral distributions as well, 

P = 1 +^^'^ P = 1 + , 
V = PP, j (fO^P = j (fe^P = 1 . (2.6) 

The expression (12.51 1 for the probability distribution represents a spurion superfield with non-zero F and D terms: 

['P]f - r ^ 0, [Vh = lei' ^ . (2.7) 

Therefore, upon averaging over the stochastic superspace, we expect to obtain a theory with softly broken F-type and D-type 
terms included. The parameter ^ then defines the supersymmetry breaking scale. 

To see explicitly how this goes consider the simplest theory with a single chiral superfield whose expansion in component 
fields is 

$(a;, e, 9) = <j){x) + V20iP{x) + 0^F{x) + iea^'Od^cjiix) + -^e^df,ip{x)a''0 - ^e^e^DcPix) . (2.8) 

V 2 4 

The Lagrangian of the model is the usual one modified by averaging over the distribution of the stochastic Grassmann coordi- 
nates, 

L={C), (2.9) 
where £ is the usual super-Lagrangian density which consists of two terms: the kinetic term, which is the D-density, 

£kin = $ + * , (2.10) 

and the superpotential term, 

W/=-$2 -^3 (2.11) 
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(plus the anti-chiral superpotential W = W'^) which is the F-density. 
The averaging of the kinetic Lagrangian ( 12.10b results in 

(C^n) = J d^eSePP^+^ = Lkin-susY + <jy* {x)(l){x) + C <l>{x)F* {x) + i(ly*{x)F{x) , (2.12) 

where Lkin-susY — J cPOdPS^^^ is the usual supersymmetric kinetic Lagrangian. Solving for the auxiliary field (ignoring 
the superpotential for the moment), 

F(x) = -r0, (2.13) 

and plugging ( 12.131 ) back into ( |2.12t we obtain that all terms proportional to ^ cancel with each other, and the remaining 
Lagrangian is invariant under the on-shell N=l supersymmetry transformations. 
Let us now compute the potential density: 

V^{W)+ h.c. = VsusY + (^(f" + ^'/>' + h.c.^ , (2.14) 

where Vsusy is invariant under N=l supersymmetry transformations. The F-term now reads: 

F = -C<t> + m*<j)* + h*<j)*^ . (2.15) 
Substituting (12.151 1 back into the total Lagrangian ( (Lkin + W)), we obtain 

/f*m 2£*h \ 

L = Lon-shcll-SUSY - ( + + ^-C- j , (2.16) 

where ion-shcii-suSY is the familiar on-shell supersymmetric Lagrangian of the interacting Wess-Zumino model. The remain- 
ing terms in ( 12.161 1 break supersymmetry softly. Thus we have a Majorana fermion ip with mass |m| and two scalar fields with 



masses y^|mp"+]^*m[ and ^\m\'^ — \^*m\. To avoid tachyonic states we must assume |^| < \m\. The supertrace STrM^ is 
zero. 

The above simple case captures some key aspects of more generic theories with stochastic supersymmetry. Namely, the 
averaged D-density, 

ikin-gauge = (Tr^+e^^^^) , (2.17) 
of a charged scalar superfield $ interacting with a gauge superfield (in the Wess-Zumino gauge), 

V{x, e, 9) = -6»a^M^(x) + ie°'0jaX'^{x) - i0j'^e°'Xa{x) + ^e^PD{x) , (2.18) 

does not produce supersymmetry breaking terms. Soft breaking terms emerge from the (gauge invariant) superpotential only.^ 

ni. STANDARD MODEL IN STOCHASTIC SUPERSPACE 

The structure of the soft breaking terms in the context of the minimal standard model with stochastic supersymmetry should 
now be clear. They all come from the F-densities. For example, averaging of the superpotential (with R-parity conservation 
being assumed), 

WsM = I^HuHa + r^QU'Hu + y'""^^QD^Hd + t^'LE'Hd (3.1) 
results in the soft breaking terms 



+ h.c. (3.2) 



where fields with tildes denote the scalar components of the corresponding chiral superfields (in a self-evident notation). Soft 
breaking masses for the gauginos A^'^ (a;) appear upon averaging of the gauge-kinetic F-densities, 



i Li 

where 



J F 



l-^TrA^A^ -hh.c. (3.3) 



are the field-strength spinor superfields for SU(3) (i ~ 3), SU(2) (i = 2) and U(l)y (i = 1) subgroups of the standard model 

(i) 

gauge group, D),, are the corresponding covariant derivatives and = x'^ + i9a^9 is the chiral coordinate. The fundamental- 
representation generators are normalized as TrT^T*' = ^6°'''. 

Summarizing, we have obtained the following tree-level soft breaking terms: 

• The bilinear Higgs soft term with = C*- 

• The trihnear scalar soft terms proportional to the Yukawa couplings, with the universal constant Aq = 2^*. 

• The universal gaugino masses, mi/2 = ^\^\- 

• The scalar soft masses are absent, TOq = 0. 

At the quantum level, the above soft breaking terms are, of course, renormalized, so the relations hold only at some energy scale 
A which joins ^ as a free parameter in the theory. 

Remarkably, the pattern of the soft breaking terms represents a particular case of the more general pattern of the so-called 
constrained minimal supersynmietric Standard Model (CMSSM). The CMSSM soft breaking terms are usually motivated by 
spontaneous supersynmietry breaking within minimal N=l supergravity.'' The model is defined through the following set of 



Gauge invariance requires the introduction of more tlian one cliiral superfield. 

In the minimal version of N= 1 supergravity one assumes the minimal flavor-blind structure for the Kahler potential, which does not follow from any symmetry 
principle. Hence, the phenomenologically desirable pattern of the soft breaking terms in the CMSSM is arguably not so well-motivated theoretically. 
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extra parameters (beyond those of the standard model): {niQ, Aq, 1111/2, tan/?, sgn/i).^ This model has been extensively 
studied in the literature and the constraints on the soft breaking parameters have been obtained from various astrophysical and 
collider experimental data. 

Obviously, our model is even more constrained than the CMSSM. Nevertheless, we shall show that it is phenomenologically 
viable in an appropriate region of A) parameter space. To do this, we solve the renormalization group (RG) equations and 
compute the sparticle spectrum at the weak scale as a function of ^, A and sgn/i. 

The RG equations are, of course, a complicated set of coupled differential equations which in general require numerical 
calculations. However, in certain regimes an analytical solution to an approximation of the one-loop RG equations is available 
We consider the regime of small/moderate tan (3 where the only important Yukawa coupling is the top- Yukawa coupling 
Ut (the other Yukawa couplings being set to zero in the RG equations). We focus on the regime Mqut < A < Mp, where 
-^^GUT — 2 • lO^^GeV is the putative grand-unification scale and Mp ~ 10^^ GcV is the Planck mass. In this regime we have 
obtained an approximate solution for tan /3,^ 



tan/3 w 5.9 (for A = A/qut) and 7.5 (for A ^ Mf 



(3.5) 



where we have taken ^ < 0. These values of tan (3 turn out to be largely insensitive to the parameter ^. For positive ^ the values 
of tan (3 are low, tan (3 ~ 2.8 for A = AfcuT and 2.4 for A = Alp. In this case the predicted mass for the lightest CP-even 
Higgs mass fails to satisfy the LEP limit The one-loop soft scalar masses for sleptons and the first two generations of squarks 
(which with good accuracy can be considered as the mass eigenstates) are: 



TO J R 



Af|) cos(2/3)/6, 



(1.70,2.09)e^ + (4Af2, 
(1.70, 2.09)^2 „(2A/2, 

(1.60, 1.97)e' - 2{M^ - Af|) cos(2/3)/3, 
« (1.59, 1.96)^2 ^ (^^2^ _ ^,^2) cos(2/3)/3, 
m?^ « (0.12, 0.15)^2 ^ cos(2/3)/2, 

(0.12, 0.15)e - (2Af 2, - Af|) cos(2/3)/2, 
s (0.038, 0.041)^2 ^ (^j2^ _ ^,^2 cos(2/3) , 



(3.6) 



where the first number in parentheses is for A = A/qut and the second for A = Mp. The third generation soft masses are 
different due to the contribution from the large top-Yukawa coupling, 

2 2 , 2 I ^ 

2 A 

2 A 
+ — , 



TO? 



Ur 



rrir 



3 

,2 



(3.7) 



with sizable mixing in the stop sector The explicit expression for A can be found in For our choice of soft parameters at 
the unification scale, A evaluated down at the A/^-scale is: /^{Mz) ~ (—3.36, — 3.92)^^. Other parameters are found to be: 
At{Mz) « (1.82, 1.96)^, Bf,{Mz) ~ (0.23, 0.19)^ (recall, ^ < 0) and fi^{Mz) ~ (1.61, 1.85)^^ _ ^,^2/2. The physical stop 



masses are: 



tl,2 2 



(to~^ + m-J =F y^(TO?^ - '^'ij^ + 4TO^(At - ncotpy 



(3.8) 



Finally, we estimate the mass of the lightest CP-even Higgs boson in the decoupling limit, niA S> Mz, of the one-loop leading 



^ It is customary to express the modulus of /j through the known standard model parameters and the bilinear soft parameter through tan f3 = (H^) / (H^^), 

using the extremum conditions for the Higgs potential. 
* The grand unified coupling is taken as aouT ~ 1/24.3. Other input parameters are taken as follows: the running top mass mt ~ 160.4 GeV (corresponding 

to the pole mass »ti^°'° Ri 170.9 GeV j^), sin^ dwi^z) ~ 0.2315, and we also assume that /i is positive. All running parameters are evaluated at the 

Z-pole, Mz ~ 91.2 GeV. 
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FIG. 1: For the choice A = A/gut, this plot shows the masses (in GeV) of the hghtest stau (marlced as 'm-ri'), the hghtest neutrahno Cmy,i') 
the lightest CP-even Higgs boson (solid line marked as 'm/i') and the lightest stop (marked as 'mti') as functions of ^ (also in GeV). To 
signify the importance of the mixing in the stop and stau sectors we have plotted also the ^ dependence of and m^i for the case when the 
corresponding mixings are ignored [dashed lines marked as 'm/j (no mixing)' and 'mri (no mixing)']. Because the lightest superpartner here 
is the stau, the A — A/gut possibility is ruled out for the minimal model. 



logarithmic approximation 



(At- fi cot (3y {hiMtM^) + 



Ml cos^ 



3 (MlMr\ 



mi Mi,) 



(3.9) 



where 



174 GeV, /i(a, h) = ^ In I 



g{a,b) 



M|) cos(2/3)/6. Ml = m?^ 



and 



Ml 



(4M, 



by — (a-b)'^ 

J + 2{M^ — Af^) cos(2/3)/3. Remarkably enough, we are very close to the so-called 
maximal mixing scenario (see, e.g., |01), where the lightest Higgs mass is enhanced significantly due to the mixing in the stop 
sector. This can be seen from Figure[Tl where we have plotted nih as a function of ^ with mixing (solid line marked as 'm/i') 
and without mixing (dashed line marked as 'nih (no mixing)'). When the mixing is neglected nih fails to satisfy the LEP bound. 

However, a similar large mixing is present in the stau sector as well, due to the universality of the trilinear soft breaking 
parameter Aq at the scale A. For the case A = Mqut this results in the lightest stau mass eigenstate being lighter than the 
lightest neutralino for ^ < —190 GeV, as can be seen from Figure[T] Note also that the region ^ > —380 GeV is excluded by 
the upper bound on the lightest stop TOj^ > 95.7 GeV (at 95% C.L.) 1^. Thus, in the absence of R-parity violation, the lightest 
stau is a stable particle, and this is excluded by observations. 

In the above calculation we have assumed /i > 0. For negative /i, we obtain a somewhat heavier stau. However, it is still the 
LSP for the region of ^'s where the bound on the mass of the lightest stop is satisfied. Finally, as has been mentioned earlier, for 
positive ^'s (and for either sign of /i) the lightest CP-even Higgs boson is too light to satisfy the LEP bound. We conclude that 
the minimal model with A = A/qut is excluded experimentally. 

As we raise A above A/qut, the stau becomes more massive relative to the lightest neutralino. Let us repeat the above 
calculations for the other extreme case, A = Alp. The results are displayed in Fig.|2] We see that for < —300 GeV, the LSP 
is now the lightest neutralino, so the phenomenological problems of the A = Mqut case are absent. 

For a given value of ^, a phenomenologically acceptable outcome is obtained for A above a certain value that lies between 
Mqut ™d Mp. Figure[3]illustrates this point. The (^, A) region above the long-dashed line is ruled out because the stau is the 
LSP. The region above the short-dashed line is also ruled out because of the experimental lower bound on the right-handed stau 
mass. The region below both the long- and short-dashed lines is broadly acceptable, with the region between the long-dashed 
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FIG. 2: For the choice A — Mp, this plot shows the masses (in GeV) of the hghtest stau (marked as 'mri'). the lightest neutralino ('m^i') 
the lightest CP-even Higgs boson (solid line marked as 'rrih) and the lightest stop (marked as 'mti') as functions of ^ (also in GeV). To 
signify the importance of the mixing in the stop and stau sectors we have plotted also the ^ dependence of mn and mri for the case when the 
corresponding mixings are ignored [dashed lines marked as 'mh (no mixing)' and 'm-ri (no mixing)']. The lightest superpartner here is the 
neutralino for ^ < —300 GeV. We have assumed > 0. 



and solid lines depicting the parameter space where the neutralino is up to 10 GeV lighter than the lightest stau. In this regime, 
co-annihilation processes in the early universe ensure the correct cosmological neutralino dark matter abundance is obtained. 
The complete superparticle spectrum for the phenomenologically- viable example ^ = —500 GeV and A = Mp is shown in 
Fig. II 

Finally, let us note that the minimal model has to be extended in order to incorporate nonzero neutrino masses. In such 
extensions, it is very plausible that the A — A/gut choice can be phenomenologically viable. For example, one can replace 
the R-parity conservation of the minimal model by the less restrictive baryon B-parity conservation condition [8l|]. Then a new 
set of interactions are allowed which render the stau LSP unstable. The very same interactions generate neutrino masses radia- 
tively. Other potentially viable models for lower values of A are the next-to-minimal supersymmetric model and models which 
include sterile neutrino superfields. Finally, one can also consider the models where ordinary F- or/and D-term spontaneous 
supersymmetry breakings are incorporated within the stochastic supersymmetry formalism. 



rV. CONCLUSION AND OUTLOOK 



In this paper we have proposed a conceptually different approach to softly broken supersymmetric theories. We have demon- 
strated that a field theory defined on stochastic superspace is equivalent to a supersymmetric theory with a specific and con- 
strained set of soft breaking terms. We have explicitly examined the case of the R-parity conserving MSSM where the Higgs 
bilinear soft parameter B^, the universal trilinear scalar coupling Aq and universal gaugino mass are given solely in terms of 
stochasticity parameter ^, and the universal scalar mass mo is zero. At the quantum level, these relations hold at a scale A which 
thus becomes a second parameter. As a result of these features we have been able to predict tan /3, as per Eq. ( 13.51 ), to express 
low energy soft masses through ^, and to derive the entire superparticle spectrum. Even though the theory is very constrained, 
it is nevertheless phenomenologically viable for A > A/gut with ^ negative and of a few lOO's of GeV in magnitude. The 
superparticle spectra predicted by stochastic superspace can be experimentally tested at the LHC. Should those results turn out 
to be consistent with a stochastic superspace explanation, then a new perspective wiU have been gained on the as yet unknown 
fundamental theory of spontaneous or dynamical supersymmetry breaking. 




FIG. 4: The complete superparticle spectrum for the phenomenologically-viable example ^ — 



-500 GeVandA = A4p. 
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